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Abstract 



We study the possible types of the nucleation of vacuum bubbles. We 
classify vacuum bubbles in de Sitter background and present some nu- 
merical solutions. The thin-wall approximation is employed to obtain 
the nucleation rate and the radius of vacuum bubbles. With careful 
analysis we confirm that Parke's formula is also applicable to the large 
true vacuum bubbles. The nucleation of the false vacuum bubble in de 
Sitter background is also evaluated. The tunneling process in the poten- 
tial with degenerate vacua is analyzed as the limiting cases of the large 
true vacuum bubble and false vacuum bubble. Next, we consider the 
pair creation of black holes in the background of bubble solutions. We 
obtain static bubble wall solutions of junction equation with black hole 
pair. The masses of created black holes are uniquely determined by the 
cosmological constant and surface tension on the wall. Finally, we obtain 
the rate of pair creation of black holes. 
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1 Introduction 



The aims of this paper are to obtain the nucleation rate and the radius of vacuum bubbles 
and to study the pair creation of black holes in the presence of a bubble wall. The processes of 
bubble nucleations and black hole productions may mimic Wheeler's spacetime foam structure 
PQ in the very early universe, which represents the spacetime no longer smooth at the Planck 
scale. Bubble nucleation [2J (3J HI [5] and Hawking-Moss type transition [6] may occur in the 
early universe. They may play an important role in selecting out our universe or inflation [TJ. 
In addition, the creation of black holes may also play an important role, as a topology changing 
process, in the spacetime foam structure. Some bubbles may have black holes. Some of them 
may be nucleated in the presence of a black hole, which is studied by Hiscock [8] and Berezin 
et al. [9], where the black hole act as an effective nucleation center for a bubble formation. 
Others may cause pair creation of black holes. On the other hand, the string theory landscape 
paradigm has many stable and metastable vacua [TO]. So the tunneling process becomes a 
remarkable event in this framework as well as eternal inflation [TT] . 

We would like to understand the mechanism how the complicated spacetime structure could 
be created and tunneling process occur. Can the theory tell us which bubble is our universe 
if we live in one? Why the particular bubble is favored? These questions are related to 
the cosmological constant problem and anthropic constraints. True vacuum bubbles can be 
nucleated within false vacuum as well as false vacuum bubbles can be nucleated within true 
vacuum. Black hole can also be created in these situation. These events cause the complicated 
structure of the early universe. Moreover, we may live in one bubble today. Our study show 
that the important thing is the cosmological constant. We will discuss that later in more detail. 
We want to classify all possible types of bubbles. It is related to the cosmological constant in 
our framework. 

In this paper, we directly compute the nucleation rate of possible vacuum bubbles using 
thin-wall approximation. We will classify the true vacuum bubble in three classes : "small" 
bubble - "large" background, "half" bubble - "half" background, and "large" bubble - "small" 
background, the terminology will be defined later. In Ref. [12], Parke evaluated the nucleation 
rate and the radius of a bubble with an arbitrary cosmological constant. In our terminology, 
the bubble in Ref. [12] corresponds to a "small" true vacuum bubble within a "large" false 
vacuum background. We carefully examine his formula to the cases of large true vacuum 
bubbles, degenerate case, and false vacuum bubbles while restricting the background as de 
Sitter. Another method to account for the initial quantum field state as the conformal vacuum 
with gravity was proposed [13J. The authors studied the Hawking-Moss transition as a limit of 
constrained instantons and Lee- Weinberg transition [14] as a false vacuum bubble nucleation. 
The vacuum tunneling with a DBI action is considered in Coleman-de Luccia type [T3] and 
Hawking-Moss type transition [16]. A new method to calculate the tunneling wave function 
that describes vacuum decay was studied by Gen and Sasaki [17] . The tunneling using junction 
condition and shell-mediated tunneling was studied [18]. The effect of the Gauss-Bonnet term 
on vacuum decay was also studied [TO]. The mechanism of the nucleation of a false vacuum 
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bubble due to nonminimal coupling was obtained [20]. The dynamics of a false vacuum bubble 
with an effective negative tension due to the coupling was also studied in Re. [21] (for a breathing 
bubble, see Ref. [22] and dynamics of a vacuum bubble, see Ref. [23J). It was proposed that the 
observed accelerated expansion of the universe is driven by the false vacuum energy of a colored 
scalar field in Ref. [21])- Recently, Marvel and Wesley studied thin- wall instanton with negative 
tension wall and its relation to Witten's bubble of nothing [25]. The negative boundary tension 
brane was analyzed in Ref. [2E]. In Ref. [21], the authors discussed the possible instanton 
solutions in Euclidean de Sitter space. They studied tunneling solution with the action having 
only Einstein-Hilbert term and cosmological term without scalar field. Brown and Weinberg 
studied thermal derivation of the Coleman-De Luccia tunneling [28J. The authors clarified the 
meaning of the bounce solution itself as well as oscillating bounce solutions [29J. Hawking- 
Moss transition was analyzed as a thermal fluctuation [50]. Recently, the relation between 
Coleman-de Luccia and Hawking-Moss transition was studied in Ref. [31] . 

Next, we will consider black hole pair creation. The vacuum of a strong field can decay 
due to the Schwinger mechanism [22] for particle production. Particle production is one of 
the decaying methods of the background field or given vacuum energy (for recent works, see 
Ref. [33]). For the black hole creation, Gross, Perry, and Yaffe [31] obtained the nucleation 
rate for black holes in a thermal bath. The nucleation of black holes in the de Sitter universe 
with a semiclassical instability was considered by Ginsparg and Perry |35j. The pair creation 
of black holes was originally obtained by Gibbons [3B] . Wu studied a single black hole creation 
at the birth of the universe using the no-boundary proposal [37] (for other frameworks, see Ref. 
[38]). In Ref. [39], the authors exhibited the new decaying process of the cosmological constant 
by thermal production of membranes. Eventually, the membrane collapses into a black hole. 
Black hole pair creation in the early universe was extensively investigated [101 HI] and in other 
frameworks [^2l H3| HU H5] . Caldwell, Chamblin, and Gibbons studied black hole pair creation 
in the presence of a domain wall [46] using the cut-and-paste procedure, where the background 
has vanishing cosmological constant and the probability was obtained. The repulsive property 
of the domain wall give rise to black hole pair creation. The calculating method in their work 
will be used in our paper. However, our solutions give rise to the background for the black hole 
pair creation with a bubble wall more naturally. 

The outline of this paper is as follows: In Sec. 2 we investigate the possible types of the 
nucleation of vacuum bubbles. We classify true vacuum bubbles in de Sitter background into 
nine types depending on whether the bubble and the background size are small, half, or large. 
We also present some numerical solutions corresponding to large bubbles with small background. 
We compute analytically the nucleation rate and the radius of false vacuum bubbles using the 
thin-wall approximation. We extend the analysis to the large true vacuum bubbles and show 
that the formula in Ref. [12] still holds in this case. Our approach is different from those in 
Refs. [2TJ STJ without scalar field. We also consider more variety of types. The relation to the 
Hawking-Moss transition will be mentioned. Furthermore, we explore the case of degenerate 
vacua in de Sitter space. The numerical solution of this tunneling was obtained in Ref. [29J. 
This case corresponds to a special transitions of the nucleation of a true vacuum bubble and 
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a false vacuum bubble. We compare this transition with the formation of a domain wall. The 
transition rate and the radius of the bubble are presented. The deformed bulk shape causes 
the contribution of bulk part in the action. In Sec. 3 we present the solution of pair creation 
of black holes in the presence of a bubble wall. Some possible types in previous section are 
adopted as the background space for pair creation of black holes. We obtain static bubble wall 
solutions of junction equation with black hole pair. The mass of created black holes is uniquely 
determined by given cosmological constant and surface tension on the wall. Finally, we obtain 
the rate of pair creation of black holes. We summarize and discuss our results in Sec. 4. 



2 The nucleation rate and the radius of vacuum bubbles 
in de Sitter space 



Let us consider the action 



S = V^jd x 

>M 



V^hd 3 x—, (1) 

dM K 



where k = 87rG, g = detg^, and the second term on the right-hand side is the boundary term 
[48J. {/($) is the scalar field potential with two non-degenerate minima with lower minima at 
$t and higher minima at <&p, R denotes the Ricci curvature of spacetime Ai, and K is the 
trace of the extrinsic curvature of the boundary dM.. Here we adopt the notations and sign 
conventions of Misner, Thorne, and Wheeler 



The bubble nucleation rate or the decay rate of background vacuum is semiclassically given 
by T/V = Ae~ B ^ h , where B is the difference between Euclidean action corresponding to bubble 
solution and that of the background and the prefactor A is discussed in Ref. [50], that with 
some gravitational corrections in Ref. |51j . We are interested in finding the coefficient B. 

We will take 0(4) symmetry for both $ and the spacetime metric g^, expecting its dominant 
contribution [52J. The general 0(4)-symmetric Euclidean metric is written by 

ds 2 = drf + p 2 (n) [d X 2 + sin 2 X (d6 2 + sin 2 9d(f) 2 )} . (2) 

The Euclidean field equations for $ and p have the form 

^// 3p', dU , . 

$" + = —, 3 

p d® 

P" = -fp($' 2 + t/), (4) 
respectively and the Hamiltonian constraint is given by 

p ' 2 -l-lf( I$' 2 - u] = 0, (5) 



3 V2 
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= ' T 
ri=o arj 



0, p v=0 = 0, and p v=Vmax = 0, (6) 



where the prime denotes the differentiation with respect to 77. We will consider only the case 
with initial de Sitter space. The boundary conditions for the bounce are 

drj 77=1) UY/ ri=ri„ 

where r] max is a finite value in Euclidean de Sitter space. The asymptotic value of $(77) is given 
by linijj-^^ $(77) ~ &f/t, where $i? is for a true vacuum bubble nucleation and $y is for a 
false one. Because of the finiteness of r] max , $(r] max ) is exponentially approaching to but not 
reaching $f/t- 

Now, we assume the thin-wall approximation scheme to evaluate B. The validity of that 
in the case of a true vacuum bubble has been examined [53]. The Euclidean action can be 
divided into three parts: B = B in + B wa u + B out . The configuration of the outside of the wall 
will not be changed before and after bubble formation. Thus, B out = 0. So we only need to 
consider the contribution of the wall and inside part. On the wall p[], B wa u = 2n 2 p 3 S , where 
S is the tension of the wall. The action has a stationary point, which determines the radius 
p of the bubble. The contribution from the inside part Bi n will be computed in the following 
subsections and will depend on whether it is true, false, or degenerate vacua. 



2.1 True vacuum bubbles 

In this section we consider possible types of true vacuum bubbles. We consider the only 
Uf > so that the exterior geometry of the bubble will remain to be de Sitter (dS) space. 
The true vacuum bubble can be classified according to the interior geometry and the size of 
the bubble. The interior geometry of a true vacuum bubble can be different depending on 
Ut = 0(flat), Ut < O(AdS), and Ut > 0(dS). The schematic diagrams of twelve possible types 
of solutions with different shapes are illustrated in Fig. [TJ The figures of the first row (1-1) - 
(1-3) in Fig. [I]- A represent flat interior geometry, those of the second row (2-1) - (2-3) represent 
AdS interior geometry, and those of the third (3-1) - (3-3) and the fourth row (4-1) - (4-3) in 
Fig. [T]-B represent dS interior geometry, respectively. 

The cases in Fig. [T]-B don't have stationary point in the action. Hence these solutions are 
not possible. The result is consistent with that of Ref. [2Tj . The authors analyzed the solution 
using the Israel matching condition [55] . 

The de Sitter background in Fig. [I]- A can be naturally called "small" for the first column, 
"half" for the second column, and "large" for the third column. As for the true vacuum bubbles 
we will call the bubble to be "small" , "half" , and "large" if the de Sitter background is called 
"large", "half", and "small", respectively. With this terminology, the three diagrams (1-1), 
(2-1), and (3-1) in the first column correspond to the small bubble or large background, those 
in the second column correspond to half sized bubble or background, and the three diagrams 
(1-3), (2-3), and (3-3) in the third column correspond to the large bubble or small background. 

The small bubble cases of (1-1), (2-1), and (3-1) were obtained in Refs. (U [12]. The cases 
of (1-2), (1-3), (3-3), (4-1), and (4-3) are considered in Refs. [27J H7], with different action from 
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ours, as a mechanism of reducing the cosmological constant [51] • There is no scalar field in their 
action and hence the junction condition was imposed on the bubble wall. In those papers, the 
magnitude of bubble wall's tension determines the type of solution. But we keep the tension 
constant in this paper. 

We now compute the contribution from inside the wall B in . The expression depends on the 
size of the bubble and the geometry. For the small bubble, it is given by 



B, 



4^ 



127T 2 



K 



dp 



P 6 Uf ~ ¥ 



(i - %p>\j T yi 2 (i - fp 2 ^) 1 / 2 



[1 - f U T p 2 f/ 2 - 1 
Ut 



(1 - iu F p 2 ) 3 / 2 - 1 



(7) 



The formula will be used for the first and second column in Fig. [TJ 

For the general expression corresponding to the third column for large bubbles and the 
fourth row in Fig. [TJ it is important to divide the integration range separately. If the dS bubble 
size is larger than half of its dS size as in the cases in the fourth row, the integration range 
of the true vacuum region in B in should be divided into two parts. Likewise, if the size of dS 
background is large than the half of its dS size, the range should be also divided into two parts. 
For example, B^ n in the case (4-3) is given by 



Bir 



An' 



pmax 



dp 



p A U T - ^ 

(1 - fp^Tp 

P 3 U T - £ 



In" / dp 

pmax 



P A U F - % 

;i - f P 2 u F y/ 2 

P 3 U F - ¥ 



(1 - ^Ut) 1 ' 2 (1 - fp 2 ^) 1 / 2 



(8) 



where p ma x(F/T) = ^3/kU f / f for dS geometry. This can be seen by the relation 



dp = ±dr) 



np 2 U 



1/2 



(9) 



where + for < rj < r] max /2, for r] = r] max /2, and - for f] max /2 < f] < 7] max . It turns out 
that the cases (4-1) - (4-3) in Fig. [TJ don't have the stationary point in the action, allowing no 
solutions. 

We now consider each column in more detail. 

In the case of the first three rows in Fig. [TJ the form of the nucleation rate and radius 
are mainly determined by the value k [20] ■ To obtain the bubble radius, the coefficient B is 
demanded to have the stationary value with respect to p variation. There are two bubbles with 
the same radius formula because Euclidean de Sitter space has the topology of S 4 . One is the 
the small bubble (the cases of (1-1), (2-1), and (3-1)) and the other is the large bubble (the 
cases of (1-3), (2-3), and (3-3)) depending on the relation among parameters. For the former, 
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'"'2 

Uf — Ut > ^r 2 -- For the latter, Uf — Ut < The existence of the above two types of bubbles 

2 

can also be seen numerically by giving the different value of a numerical parameter, k — 
while keeping the constant surface tension, = 2/3 [20]. Half bubbles are possible under the 
condition U F - U T = ^ ■ The 

cases in (1-1), (2-1), and (3-1) have the small numerical value 
of k than that of the cases in (1-3), (2-3), and (3-3). For large bubble, fj exceeds ^y 2 - while the 
particle remains near the true vacuum state for a long time in the inverted potential. This can 
be related to the fact that fj max becomes smaller as k becomes larger. Other point of view of 
varying the surface tension is in Refs. [251 [271 E] • Although the bubbles in the first and third 
columns have the same formula for the radius, they have different numerical values in general. 

The first column in Fig. [I]- A for the cases of the cases (1-1), (2-1), and (3-1) corresponds 
to the small true vacuum bubble or the large false vacuum dS background. The small bubbles 
occur when Up — Ut > 3kS 2 /A. This condition corresponds to 1 > (<^) 2 - The integral range 
need not to be divided as in Eq. (j7j) and the results for these cases were already obtained by 
Parke in Ref. [12] 

, = pI B 2B o[ {i + (£r 1 ) 2 }-D 1/2 ] (1Q) 

D ' - l ) D1/2 \ ' 

where D = [l + 2(^) 2 + (^) 4 1 , A 2 = [3/k(U f + U T )}, and A 2 = [3/k(U f - U T )\. p = 3S Q /e 
and B Q = 277r 2 S^/2e 3 are the bubble radius and the nucleation rate in the absence of gravity, 
respectively, e = Up — Ut- Here we used (1 — ^UfP 2 ^ 3 ^ 2 



[i-(fe) 2 ] 3 



3 U FP ) - ^372 • 

Case (2-1) and case (3-1) correspond to the small AdS and dS bubble, respectively. Case 
(1-1) corresponds to the small flat bubble. The radius of the bubble and the nucleation rate 
can be simplified after substituting Ut = into Eq. (flUj) as 

P= l + {p /2\r B= [l + (p /2A) 2 ] 2 ' (U) 

where A 2 = 3//ce. The results was first obtained in Ref. [1]. 

The second column corresponds to the half true vacuum bubble i.e. the half false vacuum 
background. These cases can be obtained from Parke's results by taking the limit of the radius 
of the bubble p equal to the size of the dS background, p = p max (F). This half bubble condition 
gives Up — Ut = 3kS 2 /4. This condition corresponds to 1 = (^-) 2 - Under this condition, D 
defined below Eq. (|TDJ) can be simplified as D = (S 2 + ^U T )(4/S 2 ). 

If we substitute p = p max (F) into Parke's formula, the radius of the bubble and the nucle- 
ation rate can be evaluated as follows 



2 2B 

P = , B = = r 



(12) 



Case (2-2) and case (3-2) correspond to the half AdS and dS bubble, respectively. Case 
^1-2) corresponds to the half flat bubble. The radius of a bubble and the nucleation rate can 
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Figure 1-A 



pOi) 



(4-1) PS - PS 



P0l) 



(4-2) PS - PS p(ti) 






Figure 1-B 



Figure 1: The schematic diagram for 12 possible types of true vacuum bubbles or matching with the thin- wall 
approximation. The f\ indicates the location of the wall. All the nine cases in Figure 1-A are possible solutions. 
The cases (4-1) - (4-3) in Figure 1-B don't have the stationary point in the action, allowing no solutions. 
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Figure 2: Numerical solutions for several types of true vacuum bubbles. We take k ~ 0.3849 for all cases. 
The first three figures are: (A) large flat bubble - small dS background corresponding to (1-3) in Fig. [TJ (B) 
large AdS bubble - small dS background corresponding to (2-3); (C) large dS bubble - small dS background 
corresponding to (3-3). In Fig. (D), we take e = 0.07. We can see from these figures that the radius of a bubble 
becomes larger as e becomes smaller. 



be simplified after substituting U F = into Eq. ( jl~2l as 



P 



S k 



B 



57T 



(13) 



Finally we consider the third column of Fig. [I]- A corresponding to the large true vacuum 
bubbles i.e. the small false vacuum background. For the cases (1-3), (2-3), and (3-3) in the 
third column in Fig. [H the integral range for the background should be divided into two parts, 
hence Bi n is then given by 



Bir 



4tt 2 
12tt 2 



p 3 U 



T 



3p 



%T 2 d P 



p 3 u f 



3p 



^P 2 u F yi 



p 3 u f 



3p 



;i - %p 2 u F yi 2 



dp 



1 - ^U T p 2 f' 2 - 1 (1 - f U F p 2 f' 2 - 1 (1 - f U F p 2 f/ 2 



Un 



U* 



Up 



(14) 
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This is clearly different from Eq. (J7J). 

The critical bubble size p can be obtained by minimizing B = Bi n + B wa u. The large bubbles 
occur when Up — Ut < 3kS 2 /4:. This condition corresponds to (t^-) 2 > 1. Case (2-3) and case 
(3-3) correspond to the large AdS and dS bubble, respectively. It turns out that the nucleation 
rate and the radius are given by the same form as in Eq. ffTUl) even though B in is different. Here 

we use (1 — ^Ufp 2 ) 3 / 2 = 2A ^ 3 / 2 — . Case (1-3) corresponds to the large flat bubble. The radius 
and the nucleation rate of a bubble have the same form as in Eq. ( TIT]) . 

To summarize, the results of the first column in Fig. [T] were already obtained in Ref. [12] 
and the results of the second column can be obtained from Eq. ffTUj) using the half size condition 
Up — Ut = 3kS 2 /A. The results of the third column have the same form as in Eq. (fit)]) even 
though Bi n is different. We concluded that there is only one formula by Parke for the possible 
types of the true vacuum bubble regardless of the size of the bubble or different B in . 

The true vacuum bubbles corresponding to the first three rows in Fig. [1] exist beyond 
the thin-wall approximation. As some examples, we present numerical solutions in Fig. [2] 
corresponding to (1-3), (2-3), and (3-3) in the third column of Fig. [fl-A. 



2.2 False vacuum bubbles 



We now consider the nucleation of a false vacuum bubble in de Sitter background. This was 
first studied by Lee and Weinberg [TJ|. We study four cases. The three cases will correspond 
to the third row in Fig. [1] with 77 — > r] max —77. The fourth case will be Hawking- Moss transition. 

The first case corresponds to the reflected diagram of (3-1) in Fig. [TJ The authors in Ref. 
pH] have obtained the nucleation of a false vacuum bubble within the true vacuum background 
and obtained the ratio between the decay rate of the true vacuum and that of the false vacuum. 
In this work, we assume the thin-wall approximation is valid. Using the approximation, we 
obtain the nucleation rate and the radius of a false vacuum bubble. The contribution from 
inside the wall is given by 



Bin 



S\in) - S T {m) = 4n 2 



pmax 



dp 



Air 2 



12tt 2 



dp 



p 3 U f 





Zp 



3p 



P 6 Uf - * 



3p 



(1 - fp 2 u F )v* (1 - %p 2 u T yi 2 



;i - % P 2 u F yi 2 (1 - % P 2 u T yi 2 



:i - § u T p 2 f /2 + 1 (1 - § u F p 2 f/ 2 + 1 



Un 



where we used p 2 max {F/T) 



kU 



F/T 



u h 



Thus the action becomes 



B 



12tt 2 



f U T p 2 f 12 + 1 



Un 



;i - § u F p 2 f 12 + 1 

Up 



(15) 



2n 2 p 3 S 



(16) 
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The tension of the wall is positive for the pure Einstein gravity as in the case of the true vacuum 
bubble. This case occurs when Up — Up > 3kS 2 /4. 

The formula for radius turns out to be the same as the case (2-1) in Fig. [TJ We obtain the 
coefficient B at this extremum 



B 



ggo[{l + (|^) 2 } + D 1/2 } 



(17) 



[i-(fe) 2 ] 3 



Here we used (1 — ^Upp 2 ) 3 ^ 2 = - — ' . As a result, the nucleation rate for the false vacuum 
bubble is given by changing —1 into +1 in the numerator in Eq. ([TO]) . The ratio in Ref. [H] 
can be obtained from Eq. ffTTl) and Eq. ffTUj) . 

The second case corresponds to the reflected diagram of (3-2) in Fig. [TJ This case occurs 
when Up — Up = 3kS 2 /4. The radius of the bubble has the form as in (2-2) and the nucleation 
rate is evaluated to be 



B 



2B n 



8U F 



+ D l/2 



(1* 



This result can be obtained if one substitutes p = p max (F) into Eq. f[T7j) . 

We now consider the third case corresponding to the reflected diagram of (3-3) in Fig. [TJ 

This case occurs when Up — Up < 3kS 2 /4:. The radius of the bubble has the form as in (2-1) and 

[( p0 ) 2 i] 3 

the nucleation rate has the same form as in Eq. ([T7]) . Here we use (1 — ^Upp 2 ) 3 ^ 2 = 2A ^3/ 2 — . 

We concluded that there is only one formula for the possible types of false vacuum bubbles 
regardless of the size of the bubble or different Bi n as Parke's formula. 

The fourth case corresponds to Hawking-Moss transition. If we change the upper limit of 
integral in Eq. ([16]) into p = p m ax and eliminate the contribution of the wall we can easily get 
the result obtained by Hawking and Moss. The thin-wall approximation is not considered in 
this transition. It describes the scalar field jumping simultaneously at the top of the potential 
barrier. The properties of Hawking-Moss transition were studied by many authors 



2.3 The case of degenerate vacua 

It is possible that the tunneling occurs via the potential with degenerate vacua in de Sitter 
space. The numerical solution of this tunneling was obtained in Ref. |29j. This tunneling is 
possible due to the changing role of the second term in Eq. ([3]) from damping to accelerating 
during the transition. This case corresponds to a special transitions of the nucleation of a large 
true vacuum bubble and a false vacuum bubble in dS background. The solution corresponds 
to the particle starting near at one of two vacuum states under the influence of damping term 
in one half of the inverted potential, climbing the hill under the influence of accelerating term 
in the other half of the potential, and reaching near at the other vacuum state. It will be 
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considered as the limiting case of (3-3) in Fig. [TJ This solution is the instanton solution rather 
than the bounce solution in dS background. 

Furthermore this transition is somewhat different from the formation of a domain wall. 
Domain walls can form in any model having a spontaneously broken discrete symmetry. The 
dynamical solutions with a thin domain wall was obtained by Vilenkin [H7| and by Ipser and 
Sikivie [58J. In general, one-axis in cartesian coordinates is perpendicular to the domain wall. 
Thus, z-axis = constant hypersurfaces are isometric to (2 + 1) -dimensional de Sitter space. The 
energy-momentum tensor is localized at the domain wall surface, z = 0. Thus, it behaves like 
a delta function. The (z, i)-part is a (1 + l)-dimensional Rindler space: 



dz 2 + (1 - 27iGa\z\) 2 [-dt 2 + e^ Gat (dx 2 + dy% (19) 



where a is the surface energy density. This metric has an event horizon at z = ±l/2irGa. In 
other coordinate system, the metric becomes 

ds 2 = (-^—^ drf + V 2 [-dt 2 + e^ G °\dx 2 + dy% (20) 

where Vilenkin's solution consists of two copies of the spacetime for 77 < l/27rGcr, joined 
together along the wall at 77 = ±l/27rG<7 J59J [601 ED] . It has the repulsive property providing a 
mechanism for a tunneling process, in particular a pair production of black holes. 

In our case, the solution has 0(4) symmetry both inside and outside the wall, the location 
of the wall is at rj = constant. The existence of the wall doesn't affect both the inside geometry 
and outside one because the process doesn't change the outside geometry after the transition. 
Thus, the geometry depends only on the vacuum energy. The spacetimes both inside and 
outside of the wall are a 4-dimensional spherical Rindler type in Lorentzian signature 

ds 2 = dr] 2 + p 2 (r])[-dt 2 + cosh 2 t(d8 2 + sin 2 flc/0 2 )], (21) 

where p{rf) = y^sin \J~^ri for dS, pijf) = 77 for flat, and pijf) = ^/|- sinh \J~^T] for AdS space. 
If the thin-wall approximation is assumed the transition rate has the form 



K z 



;i - § Uqp 2 ? 12 
u 



+ 2n 2 p 3 S , (22) 



where U is minima of the potential in dS space, the first term is from the contribution of bulk 
part due to the deformed, diminished, bulk shape and the second term is from the contribution 
of the wall. 

The nucleation rate and the critical radius of the wall are given by 

p= ; 2 and B = U ^ S ° (23) 
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< Ut < Uf 



Ut<0 



Figure 3: The schematic diagram of the bubble geometry with the wall. The left figure corresponds to the 
case of one half of the deformed de Sitter space and the other half of the deformed de Sitter space with same or 
different vacuum energy. The right figure corresponds to the case of one half of flat or anti-de Sitter space and 
the other half of the deformed de Sitter space. 

The location of the wall is smaller than rj max /2 because the process deforms, or diminishes, the 
shape of S" 4 . 

As for the consistency check, one can show that Parke's formula for the large true vacuum 
bubble (3-3) and Eq. (TTTT) for the reflected diagram of (3-3) become Eq. f[2"3"j) for the degenerate 
case if the condition Up = Ut = U Q is substituted into Parke's formula and Eq. ffTTl) . In other 
words, nucleation rates and radii of the large true vacuum bubble and false vacuum bubble are 
continuously connected at the degenerate case as we expected. 



3 Pair creation of black holes in the presence of a wall 

In this section we study the pair creation of black holes in the background with the wall. In 
Ref. [16] the authors studied pair creation of black holes in the presence of the domain wall using 
Euclidean junction condition with the Z 2 symmetry. The authors obtained the probability of 
the pair creation in the background with vanishing cosmological constant via the cut-and-paste 
procedure. In this work we will use vacuum bubbles obtained in Sec. [2J The Fig. [3] shows the 
schematic diagram of the bubble geometry with the wall. The geometry of outside the wall 
represented as the right half in Fig. [3] is dS, while inside geometry represented as the left half 
in Fig. [3] is AdS, flat, or dS depending on the cosmological constant. We define rj' = r\ max — r\ 
in the right half region. 

Let us consider the action 



Se= y/gEd X 
I M 



2k 2 



h E d 3 x— + a (b \/h E d 6 x, (24) 



£ K JT, 
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where the last term is a Nambu-Goto-type action on the wall. Actually, the bulk park has 
two distinct spaces with boundaries S + and respectively : the inside (or left) and outside 
(or right) of the wall. Following the work of Chamblin and Reall we can obtain the 
Israel junction conditions as well as the Einstein equations, and scalar field equation from 
the above action. From the Einstein equations, the curvature scalar is obtained as Re = 
K[V ft| I ) V a $ + 4Z7(<&)]. From the junction conditions, the trace of the extrinsic curvature is 
obtained as K = |kct. The direction of two normal vectors is taken to be pointing outward 
[60, 64, 65J. In our cases two normal vectors point in the opposite direction. After rearranging 
the terms, the action becomes 

S E = ~i \[h~E<$x - [ y/g^d 4 xU(^). (25) 



2 JT, JM 

If we think over the action Eq. ff24|) as a plausible initial setting, we have to take into account 
the contribution from the relation between the boundary term and the last term. Similar type 
of the action also appears in the domain wall production. The action is different from that of 
usual bubble nucleation [4] , which is slightly deviated from the case in flat spacetime [21 [3] . In 
the thin wall, the scalar field varies continuously between the true and the false vacuum states. 
Thus, the boundary term doesn't contribute to the action. 

We now do the following coordinate transformations suitable for 0(3) invariant configura- 
tions 

R = a/3/A sm(^/ A/3rj) sinx, T = a/3/ A tan -1 ^tan(^/ A/377) cosxj , for dS 

R = r] sin x, T = 77 cos x, for flat (26) 

R = v /37Asinh(v / A73i7)sinx, T = A/3/Atanlr 1 (ta,nh(yfA/3ri) cosx) , for AdS. 

Then the metric becomes 

ds 2 = (l ± dT 2 + fl R2) + R*dni (27) 

where the sign + in the parentheses denotes AdS space and — denotes dS space. 

The geometry of dS for this work can be divided by two parts: the inside (or left) and 
outside (or right) of the wall. Two halves of the Euclidean dS spaces described by Eq. f l27|) are 
joined at the wall. Each one has the range < R < and < T < tta/sTA. We take the 

de Sitter period f3 = 27Ta/3/A and divide by 2 into the action. 

For the right half, we use rj in Eq. (1261) instead of 77. The positions 7/ = and 77' = 
correspond to R — 0. Actually, the radii p of the inside wall and outside are same regardless of 
the vacuum energy in the cases we will consider. The radius p is related to R (the location of 
the wall without black holes), as we see from Eq. (126]) . The induced metric on the wall is given 
by 

dsl = dr 2 + r 2 (r)dn 2 2 . (28) 



13 



We consider several types of the configurations in Fig. [T] as the background space for black 
hole pair creation. In this work the general Euclidean junction condition becomes 



1 ± — r 2 — f 2 + \ 1 —r 2 — r 2 = AnGar, 

r 3 V r 3 



(29) 



where S Q = a and • denotes the differentiation with respect to the proper time measured by 
the observer moving along with the wall. The signs (— ) and (+) as a subscript of A represent 
left and right spacetime, respectively. After squaring twice, we define the effective potential for 
Euclidean junction equation to be 



V. 



eff 



1 1 

2 ~ 2 



2ttG 



a 



(A+ ± A_ 



=F 



A. 



GM 



(30) 



24nGa J ' 3 

with total energy is 0. The static bubble wall solution satisfies the following conditions 

dV ft 

V eff (r b ) = and -£^L = 0. (31) 



dr w b 

The solution can exist at r& = 3GM. The masses of created black holes are uniquely determined 
by the given cosmological constant and surface tension on the wall: 



M 



3GW3 



lixGo 



(A+±A_) 
2AnGcr 



(32) 



In this process the location of the wall with black holes is given by 
r b = 3GM or r b = 



(Mia + 



(33) 



while the location without black holes is given by 



1 



Twob 



2nGa + 



(A+±A-)V _,_ A_ 



24wGa 



(34) 



What is the effect of the black holes to the location at the wall? Before the creation, the 
location is approximately at de Sitter horizon but inside of that. After the creation, the location 
is moved toward the black hole horizon. 

The action for the production takes the form 



s E = s 



(B) 



s 



(°) 

E ' 



(35) 
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where and denote the action with black holes and without black holes, respectively. 
Note that the nucleation rate of the geometry with black holes is greater than one 
without black holes. 

(i) de Sitter with degenerate vacua 

Firstly, we consider the case of degenerate vacua in de Sitter space. The action is given by 

S { ° } = / Vhd 3 x-U Q dr d9 d<j> I drr 2 sm6 

2 J Jo Jo Jo Jo 

— ~2 7[2 Pwob~'^Q^~^ r l]ob^ (36) 

where p(r]) = ^/|- sin y~^r] for dS, p(rj) = r] for flat, and p(r]) = y^sinh y~^rj for AdS space, 

A = 8itGU , and the de Sitter period (3 = 27r-^/3/A. We divided the action by 2. Our results 
with vanishing cosmological constant are reduced as those in Ref. [IB]- We now compare the 
location of a bubble wall p with r wn u = . A 1 = in Eq. ([23]). These two values are the same 

WOb ^1+4^2(^2^2 ^ V V 



at x — 7r /2 and r = (see Eq. 

The action for the pair creation is evaluated to be 



S r = Sf-^/l- — [ T dr r d6 I** d<trr 2 smO 



2 V r 3 

•2tt 



Jo Jo 



U„ I dr I de I dj) I drr 2 sin e - S { ° ] 



JO Jo Jo 

A 3 \ 2nr hy /l - (9AG 2 M 2 )V3 
6G r 7 |l-Ar 2 | 



S C R - [ ^ 2 Gc 2 r\ + 



+ ^VL& + ^y^ r L&> (37) 

where S% represents the contribution due to the conical deficit angle at the black hole horizon, 
which has to be added in the action Se- The contribution on the conical angle was studied in 
Refs. [351 [33 ST]. 

We adopt the periodicity of the time coordinate, by the Bousso-Hawking normalization [JT] 



for Schwarzschild-de Sitter (SdS) space, is (3 — ^V^ 9 ^ m ) ^ [66] . rh denotes 

tb for the black hole horizon and rc for the cosmological horizon, respectively. In particular, 
this periodicity approaches a constant value given by (3^ = ^= in the Nariai Limit. The 
statement is related to selecting of the preferred observer having the zero acceleration. In their 
framework, the normalized surface gravity is obtained and the Bousso-Hawking temperature is 
found. There are recent studies on the thermodynamics of the SdS by Myung [67], Nariai class 
by Cho and Nam [68], analysis of SdS by Choudhury and Padmanabhan [69], other black holes 
by Park [70], and black holes in higher-order theories [71] . 
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(ii) the case (1-3) in Fig. [T] 

This case is that the left of the wall corresponds to flat space and the right corresponds to 
de Sitter space. The action turns out to be 



a 



2GM 



de 



2 77 



— SEipulk) — S E 



(o) 



A-kGM ( „ 1\ 2nr h Jl - (9A + G 2 M 2 ) 1 /3 
= h AnGar b = : T . 



SZ - norl 



A+ 3 / 27rr hV /l - (9A + G 2 M 2 )V3 \ a 



V3 



l-A + r 2 | 



2 ^ Pwob W 2 r wobi 



71 A 



12G 



|l"A + r 2 | 



sin# 



(3? 



where the action for flat space is obtained as in [16]. 
(iii) the case (2-3) in Fig. Q] 

This case is that the left of the wall corresponds to anti-de Sitter space and the right 
corresponds to de Sitter space. The action is evaluated to be 



S f 



H^ + ^ rM dT _ 

r 3 



o 



+ <7\1 



2GM A 



— -r 2 / dr. 

3 



\ P7Y p27Y 

j J d9 J d(j)r 2 sin e - S E (bulk) - s { ° ] 



1 A_ 2 ^ r h\ll + ^r 2 b 



l + A-r 2 h 



AirGarb 



1 A_ 2 \ 27ny/l - (9A + G 2 M 2 ) 
3 + 3 ^ 



1/3 



|l-A + r 2 | 



12G 



A- 



27rr ftV /l - (9A+G 2 Af 2 )V3 



ll-A+ril 



- A_ 



47^/1 + ^r 2 
1 + A_r 2 



° 2-3 ^ /A+ 3 A_ _ 4 



(39) 



where /? 



b • 3 'b . 

as m Kef. 



(iv) the case (3-3) in Fig. Q] 

This case is that the left of the wall corresponds to de Sitter space and the right corresponds 
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to de Sitter space with different vacua. The action turns out to be 
Se 



a 
4 



S C E -- \ 



2GM A_ f T - 

r z / ar_ 

r 3 U 



r 3 



o 



J J d9 J d(f)r 2 sin 9 - S E (bulk) - S { ° ] 



Se - wr\ 



( 



1 A_ 2 2irr h _y/l- (9A_G 2 M 2 )V3 



|l-A-rL 



1 A_ 2 \ 2vrr h+v /l - (9A+G 2 M 2 )V3 



ll-A+r? 



12G 



2vrr/ l+A /l - (9A + G 2 M 2 )V3 2nr h _ y/l - (9A_G 2 M 2 ) 1 /3 
— : ^ — ; h A_ 



11-A+rL 



|l"A-r 2 _ 



(40) 



In this section we have obtained pair creation rate of black holes in the four cases with different 
background types. 



4 Summary and Discussions 

In this paper we studied the possible types of vacuum bubbles and calculated the radius and 
the nucleation rate. We have obtained some numerical solutions as well as analytic computation 
using the thin- wall approximation. We consider the only Up > so that the exterior geometry of 
the bubble will remain to be de Sitter (dS) space. There are nine types of true vacuum bubbles, 
three false vacuum bubbles, and Hawking-Moss transition in which the thin-wall approximation 
is not considered. We confirmed that Parke's formula is applicable to all types of true vacuum 
bubbles even though B in is different. In addition, we obtained the single formula that applies 
also for all types of false vacuum bubbles. There are some conditions for classifying the bubbles. 
The conditions are Up — Ut > 3kS' 2 /4 for small bubbles, Up — Ut = 3kS' 2 /4 for half bubbles, 
and Up — Ut < 3kS' 2 /4 for large bubbles. The large true vacuum bubble corresponds to the 
case (3-3) and false vacuum bubble to the reflected diagram of (3-3). The coefficient B and the 
radius of the former is continuously connected to those of the latter. 

We have discussed the special case corresponding to the tunneling of degenerate vacua in de 
Sitter space. This tunneling is possible due to the changing role of the second term in Eq. (jSj) 
from damping to acceleration during the transition. This solution does not have any negative 
mode, hence corresponds to an instanton solution rather than a bounce solution in de Sitter 
background. We have compared this transition with the formation of a domain wall. We have 
obtained the transition rate and the radius of a bubble. This solution can be applied to the 
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sine-Gordan model with periodic potential. Moreover this solution can be applied to the case 
with a complex scalar field, which has global U(l) symmetry. In the potential, the transition 
between different minima, by moving along the circle of minima, requires an infinite energy 
because of the infinite volume of the space. The kinetic term in the energy is proportional 
to the volume. Therefore, a change in the field in the whole space requires an infinite energy 
and thus the transition is impossible. However, the transition via tunneling can be possible 
in de Sitter space. In addition, this transition by moving along the circle will be possible 
because of the finite volume of de Sitter space. The process of the symmetry restoration will 
be also possible via the Hawking-Moss transition. It will be interesting if this transition can 
be embedded in the case with SO (3) model as well as in the cases with gauge fields in various 
dimensions. 

We considered the pair creation of black holes in some of these vacuum bubble backgrounds. 
These solutions give rise to the background for the black hole pair creation more naturally. We 
have obtained static bubble wall solutions of junction equation with black hole pair. The masses 
of created black holes are uniquely determined by given cosmological constant and the surface 
energy density or tension on the wall. In this work, the masses of black holes are same. Finally, 
we have obtained the rate of pair creation of black holes by taking the difference between the 
nucleation rate of the geometry with black holes and that without black holes. In this 
paper we considered black holes having only the same mass in pure Einstein gravity. The case 
with different masses can be applied. Considering other theories of gravity with the modified 
tension of the wall will be also interesting. 

Our solutions, even if it has simple structure, can be used to understand the mechanism how 
the complicated spacetime structure could be created in the early universe as well as tunneling 
phenomena occur in the string landscape and eternal inflation. 

All cases we considered have at least one de Sitter space. Thus, we now need to argue 
the temperature of SdS space defined by the surface gravity. Actually, the one side of the 
background with black holes has one horizon due to a black hole because the location of the 
wall for each solution is inside de Sitter horizon. In general, SdS space has two horizons. One is 
black hole event horizon and the other the cosmological horizon, respectively. The temperature 
has also two types, Hawking temperature for event horizon and Gibbons-Hawking temperature 
for cosmological horizon in the standard normalization, respectively. Thus, there was a difficulty 
determining the periodicity of the time coordinate. In other words, there is a conical singularity 
unless both singularity may be simultaneously removed. 

Actually, the surface gravity is related to the force required at infinity to hold a test par- 
ticle hovering at the horizon. In SdS, there is no asymptotic flat region and to select the 
preferred observer with zero acceleration is difficult. Thus, one uses the standard normalization 
with two different temperatures defined by two different surface gravities. In the framework 
of the Bousso-Hawking normalization, there is one point at which the black hole attraction 
and the cosmological repulsion exactly cancel out. Thus, one can define the Bousso-Hawking 
temperature and the periodicity of the time coordinate. 

If the background becomes SdS and we can't control the singular properties an extra con- 
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tribution from the boundary term due to horizons, after black hole pair creation, will appear. 
There is no cosmological horizon in our solution because the wall is located within the horizon. 
The boundary effect due to the existence of the wall appeared in Eq. is already included 
in our solution. The contribution due to the conical angle at the black hole horizon has to be 
added in the action Se- 
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